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We consider a three-terminal system with a chiral edge channel connecting the source and drain 
terminals. Charge can tunnel between the chiral edge and a third terminal. The third terminal is 
maintained at a different temperature and voltage than the source and drain. We prove a general 
' relation for the current noises detected in the drain and third terminal. It has the same structure 

as an equilibrium fluctuation-dissipation relation with the nonlinear response dl/dV in place of the 
■ linear conductance. The result applies to a general chiral system and can be useful for detecting 

' "upstream" modes on quantum Hall edges. 
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I. INTRODUCTION 



Fluctuation-dissipation theorems (FDT)i establish a beautiful and useful connection between response functions 
^ ' and noise. They have a long history beginning with the Einstein relations and Nyquist formula and culminating in 
^ Kubo's linear response theory. The standard FDT applies in thermal equilibrium only and much attention has been 
focused on its violations in nonequilibrium conditions. It became gradually clear that the FDT forms a special case 
I , of more general fluctuation theorems valid for various classes of nonequilibrium systems^. Well-known examples are 
"j^ ' the Jarzynski equality^ and the Agarwal formulcfeS. 

', The foundations of the linear response theory and the FDT are the Gibbs distribution and causality. According to 
^ ; • the causality principle, there is a fundamental asymmetry between the past and the future since the future depends 
' on the past but the past is not influenced by future events. This imposes crucial restrictions on response to any 
^ perturbations. In this paper we address chiral systems* which possess a similar asymmetry between left and right so 
I ' that what happens on the right affects what later happens on the left but not vice versa. Obviously, this may only 
, be possible if excitations can propagate just in one direction. Such chiral transport can occur in topological states 
^ ■ of matter, a primary example being a low-temperature 2D electron gas in the conditions of the quantum Hall effect^ 
(QHE). The gas is gapped in the bulk and its low-energy physics is determined by ID chiral edge excitations. In the 
simplest QHE states all edge modes have the same chirality and hence the current flows in one direction only, e.g., 
clockwise. We show that in such chiral systems a Nyquist-type formula ^ holds for the low-frequency current noise 
. and nonlinear conductance even far from equilibrium. This far-from-equilibrium FDT is different from a more general 
Agarwal formulai^ for non-chiral systems which connects quantities that do not generally have an obvious physical 
meaning and cannot be easily extracted from experiment. 

Our results apply beyond QHE. As usual in statistical mechanics, the simplest example of a chiral system comes 
ps) [ from the physics of ideal gases. Consider a large reservoir filled with an ideal gas. A narrow tube with smooth walls 
■ ' ■ and an open end is connected to the reservoir. Molecules can leave the reservoir through the tube. The projections 
(<•-») ] of their velocities on the tube axis cannot change. Hence, they can only move from the reservoir to the open end of 
the tube and the system is chiral. Imagine now that molecules can escape through the walls of the tube with the 
probability depending on their position and velocity. A relation similar to Eq. ([T]) can then be derived for the particle 
flux through the walls and the fluctuations of the fluxes through the walls and the open end of the tube. We discuss 
that relation in the appendix. The gas example is one-dimensional. Chiral systems are also possible in 2D. Indeed, 
, topological states of matter with gapless chiral excitations on a 2D surface of a 3D system are possible (e. g., Ref. 
' and related systems). In such systems, charge can propagate in both directions along one of the coordinate axes 
but only in one direction along the second axis. Moreover, chiral models can emerge beyond conventional condensed 
matter physics. For example, statistical mechanics has been used to describe traffioi. A chiral model describes traffic 
on a network of one-way streets with no parking as long as no traffic jams form. 

Chiral edge states in QHE are of particular interest. It was proposed that non-Abelian anyons exist in QHE at 
some filling factors^, such as 5/2 and 7/2. If the prediction is true this will have major implications for fundamental 
physics and quantum information technology®. However, the nature of the 5/2 state remains an open question. 
Competing tlieories predict both Abelian and non-Abelian statistics [see Ref. ^ for a review of proposed states]. 
Some of the proposed states have chiral edges and others do not. In particular, all published proposals for Abelian 
states are chiraliS. Thus, testing chirality of QHE edges is important in this contex t "'^^'^^ and the theorem ([T|) will 
be useful for that purpose. On the other hand, it is generally believed that the edges of the Laughlin states at 
ly — l/(2p-f 1) are chiral. This expectation is supported by the chiral Luttinger liquid model^ (CLL). However, 
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CLL faces major challenges from experiment (for a review, see Ref. [l3l ). For example, it cannot explain observed 
quasiparticle transmission through an opaque barrier without bunching into electrons^ii^. Thus, it is important to 
test major assumptions of CLL. One of them is chirality. Our theorem can be used for testing that assumption. Eq. 
^ has already been verifiedi^ in the limiting cases oi T ^ V and V ^ T. 
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FIG. 1: Three-terminal setup. A quantum Hall bar is connected to source S at the voltage V. Charge tunnels into terminal C. 
The arrows represent the directions of the chiral edge modes. 

A nonequilibrium FDT can be formulated for chiral systems in various geometries. Below we focus on the simplest 
geometry illustrated in Fig. 1. We consider a quantum Hall bar connected to the source (S) and drain (D) terminals. 
The impedance between the bar and the outside world is small. Long-range Coulomb forces are screened by a gate 
(this also ensures the absence of bulk currenta^i) . Excitations propagate from the right to the left on the lower edge 
and in the opposite direction on the upper edge. The system size is much greater than the magnetic length; we 
assume that the chiral edges are far apart and do not influence each other. A third terminal C is connected to the 
lower edge through a tunneling contact. The details of the contact are unimportant and our results apply no matter 
how high or low the tunneling current It into terminal C is. A voltage bias V is applied between the source and 
C. The temperature T of the source and drain reservoirs may be different from the temperature of reservoir C. Our 
results do not depend on the latter temperature or the nature of conductor C. V and T should be much lower than 
the QHE gap (otherwise, QHE is absent and the system is not chiral). We consider the current noise in terminal C, 
Sc = J dt(A/T(t)A/T(0) + A/T(0)A/T(i)), and in the drain, Su = J dt{AI-D{t)AIuiO) + A 1^(0) A 1^(1)), where /d 
is the electric current in the drain and AI — I — {I), and derive the relation 

Su^Sc- 4t|^ + 4GT, (1) 

where G is the Hall conductance of the quantum Hall bar without a tunneling contact. This is the main result of the 
paper. 

A similar formula in a different geometry was obtained in Refs. [l8l - [20l for the exactly solvable CLL model. As seen 
below, Eq. ([1]) holds independently of the integrability of a model and does not rely on CLL. Only chirality matters. 
This point is disguised in the modeli^"— since the same solvable Hamiltonian describes a chiral system with tunneling 
between quantum Hall edges and a nonchiral quantum wire with an impurity. 

We derive Eq. ([T]) below. In the next section we give a simple heuristic derivation. Section HI contains a full 
quantum proof of the fluctuation-dissipation theorem. We summarize our results and discuss their experimental 
implications in Section IV. Appendix contains a derivation of Eq. ([ij in an ideal gas model. 



II. HEURISTIC DERIVATION 



Equation ([T]) does not contain the Planck constant and so we first give its heuristic classical derivation. The current 
Id = /l — is composed of the current /l , entering the drain along the lower edge, and the current lu on the upper 
edge. These currents are uncorrelated and hence the noise in the drain is the sum of the noises of I\j and /l: 

'5'd = '5'u + 'S'l- (2) 

The noise on the upper edge is the same as in the symmetric situation without tunneling into C. In the latter case 
the noise Su is given by the Nyquist formula. Thus, S\j is one half of the equilibrium Nyquist noise, = 2GT. In 
order to evaluate we note that in a steady state there is no charge accumulation on the lower edge and hence the 
low-frequency component of the current, absorbed by the drain, Jl = /s — /t, where Ig is the low-frequency part of 
the current, emitted from the source. Thus, 



Sh = Sc + Ss — 2SsT, 



(3) 
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where the cross-noise Sst = J dt{AlT{t)AIs{0) + A/s(0)A/T(i)) and the noise of the emitted current equals one half 
of the Nyquist noise because of chirality, 

Ss=Sv = 2GT. (4) 

We are left with the calculation of the cross-noise. The tunneling current depends on the average emitted current GV 
and its fluctuations I^j- We assume that the central part of the lower edge has a relaxation time t. It is convenient 
to separate the fluctuations of 7s into fast, />, and slow, parts. contains only frequencies below 1/t. An 
instantaneous value of the tunneling current /t(0 depends on the emitted current within the time interval r. does 
not exhibit time-dependence within such time interval. Hence, it enters the expression for the tunneling current in 
the combination GV + only. It — {I{GV + /^)), where the brackets denote the average with respect to the 
fluctuations of /g. According to the Nyquist formula for the emitted current, its harmonics with different frequencies 
have zero correlation functions: {luji-u)') ~ 5{uj—uj'). For the sake of the heuristic argument we will assume a Gaussian 
distribution of Is and hence independence of its high and low frequency fluctuations (no such assumptions are needed 
in a general proof). After averaging with respect to the fast fluctuations of Is we can write It = {J{GV + I'^)), 
where J is obtained by averaging over corresponds to a narrow frequency window and can be neglected in 

comparison with GV, i.e.. It = J{GV). For the calculation of the cross-noise we expand J{GV + to the first 
order in and obtain 

Sst = {It,uiI-lj + It.-^Ilo) = qqy {^^^-^ + I-i^Iuj) = 2T-^, (5) 

where we used the Nyquist formula for the fiuctuations lu: oi Is- A combination of Eqs. I^j^ yields the desired result 
([1]) for nonequilibrium steady states. Note that only low-frequency particle current conserves along the edge in our 
screened conductor and hence it is plausible to expect that It depends only on harmonics luj with hui < fiwcutoff ^ T. 
In such case the assumption of the Gaussian distribution and independence for is not needed since the same 
results can be obtained from the lowest order expansion in 

The above argument easily generalizes to other geometries with many terminals and/or tunneling between QHE 
edges. 

One can estimate noises and Sc- We assume that V ^ T ^ lOOmK. These are realistic parameters for noise 
experiments. The injected current Is ~ e^V/h since G ^ e'^/h. We assume that the transmission probability to 
conductor C is of the order of 1/2. Thus, It ^ Is- Then the second and third terms on the right hand side of Eq. 
(1) are of the order of e'^[kT/h] ^ e'^[eV/h]- The noise Sc can be estimated from the fact that its physical meaning 
corresponds to the ratio of the fiuctuation (AQ^) of the transmitted charge to the time At over which the charge was 
transmitted to C. The probability of a single tunneling event during the time interval At ~ h/eV is of the order of 
1/2. Thus Sc ^ I At ^ e^[eV]/h and has the same order of magnitude as the second and third terms on the right 
hand side of Eq. (1). Finally, So has the same order of magnitude. This corresponds to the noises of the order of 
I H z . This is a very small number but even lower noises are measured in the state of the art experiments in 

the field. 

Counter- propagating edge modes break Eq. ([T]). Imagine that an "upstream" neutral modes propagates in the 
direction, opposite to that of the charged mode. Each tunneling event into C creates a neutral excitation that brings 
energy ~ V back to the source and heats it. This increases the noise, generated by the source, and raises the effective 
temperature above T in Eq. ([T]). We would like to emphasize that the bulk of the source remains at the temperature 
T since the heat capacity of the source is large. However, the bulk plays relatively little role in noise generation in 
our system due to a relatively low resistance of a massive bulk conductor. 



III. PROOF OF THE NONEQUILIBRIUM FDT 

We now turn to a full quantum derivation of Eq. ([1]). To simplify notations we will omit the Boltzmann constant 
and h from the equations below. As is clear from the above classical argument, the drain potential has no effect on 
the noises. It will be convenient to assume that the source and drain potentials are equal. It will also be convenient 
to assume that at the initial moment of time t = —oo there was no tunneling or other interactions between conductor 
C and the QHE subsystem containing the quantum Hall bar and the source and drain reservoirs. Thus, the system 
is initially described by the Hamiltonian i7o — He + Hu, where He and H^i denote the Hamiltonians of the two 
subsystems: conductor C and the QHE subsystem respectively. At later times the Hamiltonian includes an interaction 
term, H = Hq + II\{t)- One of the effects of the interaction is charge tunneling between the QHE subsystem and 
conductor C. We assume that the interaction Hi becomes time-independent well before the moment of time t = and 
the system is in its nonequilibrium steady state at f = 0. The steady state depends on the voltage bias V , temperature 
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T and the temperature Tq of conductor C. All these energy scales must be lower than the QHE gap. Otherwise, the 
system is unlikely to allow a chiral description. We do not assume any special relation between the energy scales. 
li T — Tc and V — Q then the system is in equilibrium. If T — Tc ^ T and V ^ T then the system is close to 
equilibrium. Otherwise, the system is far from equilibrium. Our main result ([T]) applies in all those cases. 



A. Chiral systems 



A general definition of a chiral system is the following: Consider a system whose Hamiltonian has a time-dependent 
contribution Ht = J^^dxh{x,t), where the integration extends to the left of point y. In a chiral system, local 
observables to the right of point y do not depend on the form of h{x,t) for any initial conditions. 

In what follows, we will not need the most general definition. Instead, we will focus on one particular observable, 
current /g, emitted from the source to the lower edge. We are only interested in the low- frequency regime. In that 
limit, the precise choice of the point, where the current /g is measured, is unimportant due to the charge conservation. 
The same low- frequency current flows in all points of the lower edge between the source and conductor C. Let us 
select a point A in the gapped region in the bulk of the QHE liquid and a point B below the QHE bar (Fig. 2). The 
same current Jg must flow through any line, connecting A and B. This remains true, even if the line does not cross the 
lower edge and instead goes through the source (and the boundary between the source and the QHE liquid). Thus, 
/g can be defined both in terms of the edge and source physics. 
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FIG. 2: The same low-frequency current Is flows through both dashed lines. 

The chirality assumption means that the average emitted current /g [t) does not depend on the presence of conductor 
C for any initial conditions. In other words, the expression Trp(< = — oo)/g(t), where p is the initial density matrix, 
is the same when the Heisenberg operators 



h{t) = [Texp(-z / V{t)dt)]-^Is{-^)[TeM~i f V{t)dt)], (6) 

J —OO J —OO 

where /g(— oo) is a Schrddinger operator, are defined in terms of the Hamiltonians V = Hq and V{t) = HQ+Hi{t) = H . 
The latter can be true for a general p(— oo) only if the Heisenberg operators /g(t) are the same in the presence and 
absence of conductor C. Note that this property is satisfied in the chiral Luttinger liquid model. To the right of 
conductor C, the CLL action assumes the form L — mh/ {ATT)[dt4)dx4i — w(9a;(/))^], where v is the velocity of the edge 
excitations, the charge density q — e9a;(/)/(27r) and the current operator /g — vq. From the equation of motion, 
dtdx4> — vd^cf) = 0, we find that the electric current on the right of conductor C, /g = /g(i -I- x/v), depends only on 
the initial conditions on the right and is not affected by the form of Hi{t). The same property is satisfied in any 
other chiral conformal theory and in many other situations. For example, the chirality property of the operator /g 
survives, if any changes are introduced into the above CLL action to the left of the point, where /g is measured. The 
chirality assumption also holds for QHE edges with several modes of the same chirality but breaks down, generally, 
if counter-propagating modes are present, as e.g., in the anti-Pfaffian stat o^^'^^ proposed at = 5/2. 

Edge reconstruction^^ may result in "net chiral" edges that are not chiral. For example, a pair of counter-propagating 
integer QHE modes can emerge on a = 1/3 edge. In general, this breaks chirality. However, in practice, disorder is 
likely to localize such mode pairs and restore chirality. 
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B. Initial density matrix and Heisenberg current operator 

At the time t — —oo there is no interaction between conductor C and the QHE subsystem that includes the 
source, drain and 2D electron gas. Hence, the initial density matrix p{—oo) — pnPc factorizes into a product of 
the initial density matrix pc of conductor C and the initial density matrix pn of the QHE subsystem. Each of 
them corresponds to an independent Gibbs distribution determined by an appropriate reservoir. At later times the 
subsystems interact, the steady state depends on both reservoirs, and the factorization property no longer holds in the 
Schrodinger representation. Thus, it will be convenient for us to perform calculations in terms of the initial density 
matrix because of its simpler structure^!. This means that we will use the Heisenberg formalism so that all time 
dependence is placed into the operators of observables. The chirality property will allow us to extract considerable 
information about the matrix elements of the Heisenberg operator Is {t) and prove Eq. ([T]) . We would like to emphasize 
that the Hamiltonian has a time-dependent piece Hi(t) and this piece enters the definition of all Heisenberg operators. 
The presence of that piece is crucial for the difference between the density matrices in the Heisenberg and Schrodinger 
representations. We will omit the time argument in pjj and pc- It will be always understood that these are initial 
density matrices at t = — oo. In all calculations below, p is also taken at t ^ — oo. Certainly, in the Heisenberg 
representation, the density matrix does not depend on time. 

Our approach resembles the Keldysh formalism, where all correlation functions are also expressed in terms of the 
initial density matrix. In the Keldysh technique, if the interaction is adiabatically turned on the initial density matrix 
describes free particles and hence factorizes into a product of single-particle density matrices. A difference from 
our approach consists in the application of the interaction representation in the Keldysh perturbation theory. The 
average of any properly time-ordered product of creation and annihilation operators is known exactly in the interaction 
representation. This allows development of a diagrammatic technique. We use the Heisenberg representation instead 
and rely on special properties of the matrix elements of the operator /g in the basis, in which the initial density matrix 
is diagonal. 

Any density matrix is Hermitian and can be diagonalized. Hence, 

Pii,c = ^PB.,Cn\nK,c){nH,c\ (7) 

and 

p(-oo) = ^p„|n)(n|, (8) 

where 

Pn = Pmi'PCn" (9) 

and 

= (10) 

where the states |n^) and jn'^) are selected from the Hilbert spaces of the QHE system and conductor C respectively. 
Ph is a Gibbs distribution, pHn ^ exp(— E'n/T — \e\VNn/T), where Nn is the number of electrons in the QHE 
subsystem, T the temperature of the source and drain reservoirs and En are the eigenenergies of the eigenstates 
Itih) of the quantum Hall subsystem before the tunneling contact was turned on, i.e., jrin) are eigenstates of the 
Hamiltonian Hu with particle numbers 7V„ . We do not make assumptions about pc . Our proof applies as long as the 
initial density matrix p(— oo) factorizes and the initial density matrix of the QHE subsystem pn is given by the Gibbs 
distribution. In practice, the initial density matrix pc is also likely to be a Gibbs distribution. To avoid a possibility 
of confusion, we emphasize that all states in the bases jrifj) and \nc) are time- independent. Thus, they are no longer 
eigenstates of the time-dependent Hamiltonian after the interaction Hi has been turned on. 

If the interaction Hi{t) is never turned on then Is{t) acts in the Hilbert space of the QHE subsystem and hence 
its nonzero matrix elements are always diagonal in the basis of |nc). The chirality property means that the same 
restriction applies to nonzero matrix elements of /s('^) even after the interaction Hi{t) has been turned on since Is(t) 
must be the same in the presence and absence of the interaction Hi{t). The emitted current operator commutes with 
the number N of the particles in the quantum Hall subsystem since it describes particle transfer between the source 
and the edge. Thus, in the absence of the interaction Hi{t), it has nonzero matrix elements only between states jrifj) 



6 



with the same Nn- Again, the chirahty property means that the same restriction on nonzero matrix elements applies 
even after the interaction has been turned on. Before the interaction between the quantum HaU bar and subsystem C 
has been turned on, it is easy to write the time-dependence for matrix elements of any operator acting in the Hilbert 
space of the QHE subsystem: {n\Oii{ti)\m) — exp(i[i?„ — Em\{ti — t2)){n\0'ii{t2)\m) . The same relation would apply 
at all times, if the interaction Hi(t) were never turned on. The chirality property means that the emitted current 
operator Is{t) exhibits exactly the same time-dependence at any times, if the interaction H\{t) is turned on and if it 
is not. This applies both before and after the tunneling between two subsystems has been turned on. Setting i2 = 
in the above relation, we obtain 

{n\h{t)\m) = cxp(^[S„ - i;,„]i)(n|/s(0)|m). (11) 
C. Voltage bias 

A standard way to include voltage bias in mesoscopic systems is based on the Landauer-Biittiker formalism: One 
assumes that the tunneling term is initially absent and then turned on and that the lower edge is initially at equilibrium 
with the reservoir with the chemical potential V . We will use a mixed Kubo-Landauer formalism to determine the 
response of /t to a small change 5V of the voltage bias. It will allow us to reduce the problem of nonlinear response 
to V to the linear response to 5V . In the mixed Kubo-Landauer language, an additional electromotive force 5V is 
generated by an infinitesimal time-dependent vector potential 5 A.. 
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FIG. 3: (color online). Illustration of the bias voltage. 5 A and 5E are applied in the region between two solid vertical lines. 
In the example in the figure the region with (5E crosses both the source (shaded) and the gapped QHE region (white). 5(f> is 
constant on the vertical dashed line. = in point Q and 54> = 5V in point P. 

We assume that different contacts are connected with infinite reservoirs at different electrochemical potentials. Their 
diff'erence determines the voltage bias V: the source electrochemical potential is V and the potential of conductor C 
is 0. A small change of the bias 5V can be introduced with an electric field described by a time-dependent vector 
potential, (5E — —l/cdSA/dt. The electric field is applied in a finite part of the source terminal (Fig. 3) and 
cannot affect chemical potentials of the infinite reservoirs. The chemical potentials determine electric potentials of the 
reservoirs because of charge neutrality. The distribution of charges certainly changes in the middle of the conductor 
in the presence of i5A, so the electrostatic potential </> also changes. However, the potential difference between the 
reservoirs does not. The magnetic field must be time-independent inside the sample as required by the restrictions on 
e.m.f. sources in the circuit theory. In other words, (5B = curWA — inside the conductor. Hence, the integral of SE 
does not depend on the choice of a path inside the conductor at fixed positions of its ends. If a path PQ begins in the 
infinite source reservoir and ends on the opposite side from the region with the field (5E (Fig. 3) then Jp drJE = dV. 
Hence, SA — ct x grad6(p, where 6(j) — SV in the source reservoir and 6(1) — far on the left in the quantum Hall 
region (Fig. 3). As a consequence, the vector potential 6 A can be gauged out inside the conductor at the expense 
of changing the electrostatic potential </>—>■ ^ -I- 6(1). This means a change of SV in the electrochemical potential of 
the source reservoir and no change in the potential of conductor C. Thus, one can see that the Kubo formalism is 
equivalent to the Landauer-Biittiker approach in the presence of infinite electrically neutral reservoirs. 

SA generates a correction to the Hamiltonian: SH — — J cfirSAj/c, where j is the current density. Consider an 
arbitrary surface of constant S(j) in the region with nonzero (5E (Fig. 3). Similar to the discussion of Fig. 2, in the 
low- frequency limit, the total current through any such surface is the same and equals the total current through the 
source 



I = Iv-Is, 



(12) 
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the signs in front of /u,s reflecting our conventions about current directions, Fig. 1. This aUows rewriting 

SH = ~ISA/c, (13) 

where d6A/dt = c6V. The same approach can be used to describe small changes of the drain potential but they are 
irrelevant for our purposes. 

In what follows it will be convenient to consider the case of SA oscillating with a low frequency w, SA = cSVsin ujt/uj. 

D. Main argument 

We now give a full quantum derivation of Eq. ([T]). The arguments leading to Eqs. ([SH]) do not change compared 
to Section II and we concentrate on Eq. ([5]). The cross- noise can be expressed as 



2S'sT = J dt{lTiO)Is{t) + h.c.){eyi-p{iujt) + ex.p{-iujt)) = 

dt^(exp(iwt) +cxp(-iwt))[(m|/T(0)|n)(n|/s(i)|m)/9™(-oo) + (n|/g(t)|m)(TO|/T(0)|n)p„(-cx))], (14) 



where a low frequency cj < 1/r (r is the relaxation time, Section II), /?(— oo) is the initial density matrix and /t.s are 
Heisenberg operators [see Eq. ([6])] . As usual, introducing a small nonzero frequency allowed us to write the expression 
in terms of /s,t and not A/g t = Is,t — (-^S.t)- Eq. (ITil) gives the noise at t = 0, when the system is in a steady 
state. As discussed in Section III.B it is convenient to use the Heisenberg representation in which the density matrix 
is not the steady state density matrix p{t = 0) but the initial p{t = — oo) since Is{t) exhibits remarkable properties in 
such representation. Inserting the time dependence (jlip of the matrix elements {n\Ig{t)\m), one finds 

2SsT = 2^^[p„,(-oo) +p„(-oo)](m|/T(0)|n)(n|/s(0)|m)[(5(£;„ w) +5(i;,„ -i;„ +w)]. (15) 

ran 

Next, we need to compute Rr = dlr/dV. As discussed in Section III.C, this is a linear response problem with 
respect to 6V. Similar to Ref. l25l . Rt is given by the same Kubo formula as in equilibrium. Indeed, 



(/T(i-O)) =Tr[p(-^)5A(-^,0)/45A(0,-w)], (16) 



where is a Schrodinger operator, SA{t2,ti) the evolution operator, S'a(0,— oo) ~ Texp(— i j'^^ HA{t)dt), the 
Hamiltonian HA{t) — H — ISA/c and H — Hq + Hn + Hi{t). The expansion to the first order in 5 A yields 

RTx5V = i I dtTT[p{-oo){Si-oo,t)6H'Sit,0)IiS{0,-oo)^ S{-oo,0)IiS{0,t)SH'S{t,~oo)}], (17) 

where SH^ is the Schrodinger operator and S{b, a) = Texp(— i J^^ dt[Hc + iJn + Hi{t)\). Substituting in the 
above equation we see that the response of /t to 5V expresses as the sum of the responses of /t to the perturbations 
Is6A/c and and —luSA/c. The latter response is zero since the edges are far apart and perturbations on the upper 
edge have no effect on the lower edge. With this in mind, we rewrite the nonlinear response to V in the form 

Rt - dlr/dV = t lim T dtY^ ^^^^^^^^^^^^-^^^^^t^ 

(18) 

In the above equation we absorbed evolution operators into the Heisenberg current operators. 

It is convenient to combine the above response with the response Rs of /§ to the perturbation SVI^ sin(a;t) /uj in 
the Hamiltonian. Certainly, that response is zero because of chirality. Indeed, we consider a perturbation, acting on 
the left of the point, where Is is measured. We get an expression of the same structure as above with the indices S 
and T exchanged. In a steady state we expect that (/s(0)/T(i)) = (/s(— t)/T(0)). This allows rewriting Rs in the 
form 



Rs - ^ lim / °°dty] ^^P^''^^) ^^^P( '^^^ [(n|/s(^)|m)(m|/T(0)|n)p„(-oo) - (m|/T(0)|n)(n|/s(i)l™)P™(-oo)] 



(19) 
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We next compute Rt = Rt + Rs'- 



Rt = 2tt lim > S{E.n - S„ + w) 



mn 



Pn(-Oo) - p,n{-Oo) 

2w 



[(n|/T(0)|m)(m|/s(0)|n) + (m|/T(0)|n)(n|/s(0)|m}], (20) 



where we used the time-dependence (jlip . The above equation contains the initial density matrix at time t = —oo 
and the Heisenberg current operators ^ at time t = 0. Finally we apply the results of Section III.B for Is{t) and 
p(— oo). We notice that nonzero matrix elements (m|/s(0)|7i) correspond to Nn — Nm and |ric) = \'mc)- Hence, in 
the limit of low frequencies in Eq. ((20|) . [p„(— oo) — pm(~<X))]/ijj = — pcndpYin/ dEn — p„(— oo)/T, where we used 
the factorization property ^ which is only valid for the initial density matrix. Comparison of Eqs. (jl5|) and (j20[) at 
small uj establishes Eq. ([5]). 

The above calculation relies on the structure of the initial density matrix p(— oo). This does not mean that the 
steady state depends on minor details of the initial state. Only the temperatures and chemical potentials of the large 
reservoirs are important. Those temperatures and potentials remain the same in the initial and steady state. If, on 
the other hand, one of the reservoirs is not large then the steady state does not depend on the initial density matrix of 
that reservoir. This can be easily seen from Eq. ([T]) in the limit of a small reservoir attached to conductor C. Indeed, 
in that case, 1^ — Sq — Q in a. steady state since conductor C cannot accumulate charge. Thus, Eq. ^ reduces to 
S'd = 4:GT. This is a usual Nyquist formula, valid for a system in thermal equilibrium at the temperature T and 
a uniform chemical potential. Obviously, the steady state is indeed an equilibrium state with the temperature T, if 
conductor C is attached to a finite reservoir. In this example, the final state does not depend on the initial density 
matrix pc- 

A similar argument does not work for a finite source reservoir. Indeed, the derivation of Eq. ([1]) relies on the 
assumption that /y is uncorrelated with /g . If the source reservoir is not large then the assumption is violated in the 
steady state and lu — Is instead. 



The focus of the preceding section is on QHE, but similar non-equilibrium FDT apply in many other systems. The 
simplest example of a chiral system, based on an ideal gas, is considered in the appendix. Our results can also be 
generalized beyond ID, for example, for the surface transport in a 3D stack of QHE systems. 

The geometry of Fig. 1 allows only electron tunneling to conductor C. FDT's, similar to ([1]), can also be derived 
in other geometries, where fractionally charged anyons tunnel: One can consider tunneling between two edges of the 
same QHE liquid. 

Eq. ([T]) does not contain the temperature Tc of conductor C. This, certainly, does not mean that the properties of 
the system do not depend on it. The current It and the noises S'd and Sc are all affected by the temperature of C. 
The general relation ([T]), however, remains the same. We would like to emphasize that our derivation does not contain 
any assumptions about the character of the dependence of It and Sc on the temperature and voltage. An interesting 
situation is possible, if conductor C is chiral and Tc ^ T. One can then derive two equations of the structure ([T]) 
with two different temperatures in them. 



FIG. 4: A possible experimental setup. Charge carriers, emitted from the source, can either tunnel through the constriction Q 
and continue towards the drain or are absorbed by the Ohmic contact C. 



IV. DISCUSSION 



— D 




Our main result, Eq. ([TJ, applies in chiral systems and can be used for an experimental test of chirality. A convenient 
measurement setup is illustrated in Fig. 4. Several mechanisms break chirality and can lead to the violation of Eq. 
dl]). One mechanism involves long range forces in the 2D electron gas. Our discussion assumed that a gate screens 



9 



long-range Coulomb interaction. This allowed us to assume that the tunneling Hamiltonian Hi does not depend 
on the voltage bias and the bias manifests itself in the 2D electron gas only through the chemical potential of the 
lower edge. Without screening, Hi may depend explicitly on the voltage and this must be taken into account at the 
calculation of 9/t / dV. Strong interaction of edge modes with non-chiral bulk modes may also break chirality. 

The most interesting mechanism of chirality breaking involves "upstream" modes^^'^^, Fig. 5. In the simplest 
example, two charged modes carry charge in the opposite directions. Let us imagine that the two chiral channels do 
not interact and all charge tunneling into C occurs due to particles, populating the upstream channel, directed from 
D to S. Then the noise in the drain is the same as in the absence of C, in contradiction with Eq. ([1]). 















Iv 
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II ^ ^ Is 
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FIG. 5: A non-chiral system. The solid line along the lower edge illustrates the "downstream mode", propagating from the 
source to the drain. The dashed line shows a counter-propagating "upstream" mode. 

This discussion neglects a possible heating effect. To illustrate it, let us assume that the upstream mode is neutral. 
It cannot carry charge but carries energy. In general, the tunneling operator into C includes a product of operators 
creating charged and neutral excitations on the edge. A neutral excitation of the energy ~ V travels to the source 
and heats it. This affects noise, generated by the source, and leads to the violation of Eq. ([1]). The details of the 
interaction of a neutral quasiparticle and the source are poorly understood theoretically. The experiment suggests 
that the heating effect will be strong^^. Thus, large deviations from Eq. (IT|) can be expected in the presence of 
"upstream" neutral modes. 

Our only assumption about V, T and Tc was that they are much lower than the QHE gap. Otherwise, a chiral 
description is unlikely to apply. If the system is chiral we make no assumptions about the relation between V and T. 
Nevertheless, our main focus was on the regime with V ^ T. Eq. ([1]) greatly simplifies and becomes less interesting 
in the opposite limits V and T V . In the former case, let us set T to zero. Then Eq. ^ reduces to Sy) ^ Sq- 
This relation reflects noiseless character of the emitted current. In the opposite limit, let us assume that V = Q and 
T = Tc- Then the equilibrium FDT applies. dlT/dV is now linear response. Hence, Sc — ^Tdl^ / dV . Finally, Eq. 
(H)) reduces to Sn — AGT. This simple relation reflects the fact that the lower edge is in thermal equilibrium on both 
sides of the contact with C. 

In conclusion, we established a non-equilibrium FDT ([TJ for chiral systems, both close {V <^ T) and far {V > T) 
from equilibrium. The result does not apply to non-chiral conductors and can be helpful in the search of counter- 
propagating modes on quantum Hall edge a^^i^^ . 
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Appendix A: Ideal gas model 

In this appendix we address a non-equilibrium FDT for an ideal gas system, briefly discussed in the introduction. 

We consider a large reservoir filled with an ideal gas of non-interacting molecules at the temperature T and chemical 
potential ^. Molecules can leave the reservoir through a narrow tube with smooth walls (Fig. 6). Collisions with the 
tube surface are elastic and do not change the velocity projection on the tube axis. Thus, molecules only move from 
the reservoir to the open end of the tube and the system is chiral. Imagine now that molecules can escape through a 
hole in the wall of the tube. We derive a relation similar to ([l}: 

5D = 5c-4r^ + 5s, (Al) 
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where It is the particle current through the hole in the tube wall, and Su = J ^^(A/d (<)A/d(0) + A/D(0)A/D(i)), 
S's = / dt(A/s(t)A/s(0) + A/s(0)A/s(0), and = / dt(A/c(t)A/c(0) + A/c(0)A/c(0) are respectively the particle 
current noises at the open end of the tube, at the opposite end of the tube, and at the hole (Fig. 6). The noise Ss 
can be determined from the measurement of Sd in the geometry without a hole in the tube wall. 



\Sc 
It 



t % 



FIG. 6: Ideal gas in a reservoir with a tube. 

The simplest proof of Eq. (|A1[) is a direct calculation along the lines of Ref. [13. The calculation is especially simple 
in the case of an ideal classical gas which should be understood as a Fermi gas with a high negative chemical potential 
in order to use the above reference. Quantum Fermi- and Bose-gases are also easy to consider. The current and noise 
can be expressed as sums of contributions from different small energy intervals. Let / = l/[exp({iJ — n}/T) + 1] be 
the Fermi distribution function for a particular energy and Tg the transmission coefficient through the tube wall for 
that energy {Te may depend on the channel number, if there are many channels). According to Ref. '2?, in a Fermi 
gas, the contribution from a corresponding energy window to the current, tunneling through the walls, is proportional 
to TeI, the contribution to Ss is determined by 2/(1 - /), Sc by 2r£;/(l -Tg/) and 5d by 2(1 -r£;)/(l - (1 -Tg)/). 
A combination of these contributions gives the desired theorem (jAl[) . 

One can also generalize our QHE proof. This approach, certainly, is harder than a direct calculation. The situation 
simplifies for a degenerate Fermi gas whose particles can tunnel outside the tube only for energies, close to the Fermi 
level. Such gas can be mapped onto a model of charged particles whose mutual interaction is completely screened 
by the gate. The electric current and noise of such charged particles equal their mass current and noise up to a 
trivial coefficient. The connection of the chemical potential and voltage bias is obvious. Such model would describe 
left-moving electrons in a quantum wire in the language of the Landauer-Biittiker formalism. Its low-energy effective 
Hamiltonian is related to the integer QHE edge physics. Tunneling through the tube walls plays exactly the same 
role as the tunneling into conductor C in the QHE setting. The only important difference from a QHE setting, Fig. 
1, is the absence of the upper edge. Thus, the derivation from the paper can be repeated with only one modification: 
Su should be set to zero in Eq. (2). A small modification involves then Eq. (4): now Ss simply equals the noise at 
the open end of the tube in the absence of the hole in its side. That quantity must be substituted instead of AGT in 
Eq. (1). Nothing else changes in that equation. 
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